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Abstract. We present several algorithms for computing the visibility polygon of a simple 
polygon V from a viewpoint inside the polygon, when the polygon resides in read-only 
memory and only few working variables can be used. The first algorithm uses a constant 
number of variables, and outputs the vertices of the visibility polygon in 0(nr) time, where 
f denotes the number of reflex vertices of V that are part of the output. The next two 
algorithms use O(logr) variables, and output the visibility polygon in 0{n\ogr) randomized 
expected time or 0(n log 2 r) deterministic time, where r is the number of reflex vertices of 
V. 

1 Introduction 

The visibility polygon of a simple polygon V from a viewpoint q is the set of all points of V that can 
be seen from g, where two points p and q can see each other whenever the segment pq is contained 
in V . The visibility polygon is a fundamental concept in computational geometry, and one of the 
first problems studied in planar visibility. The first correct and optimal algorithm for computing 
the visibility polygon from a point was found by Joe and Simpson [11 . It computes the visibility 
polygon from a point in linear time and space. We refer the reader to the survey of O'Rourke |14) 
and the book of Gosh [5] for an extensive survey of this problem and other interesting variations. 

In this paper we look for an algorithm that computes the visibility polygon of a given point 
and uses few variables. This kind of algorithms not only provides an interesting trade-off between 
running time and memory needed, but also is very useful in portable devices where important 
hardware constraints are present (such as the ones found in digital cameras or portable phones). 

A significant amount of research has focused on the design of algorithms that use few variables, 
some of them even dating from the 80s 12 . Although many models exist, most of the research 
considers that the input is in some kind of read-only data structure. In addition to the input values 
we are allowed to use few additional variables (typically a logarithmic amount). 

One of the most studied problems in this setting is that of selection. For any constant e € (0, 1), 
Munro and Raman (T3] give an algorithm that runs in 0(n 1+e ) time and uses 0(l/e) variables. 
They also give a randomized algorithm whose expected running time is O(nloglogn) under some 
assumptions on the input. Frederickson [8] improved this result with an algorithm whose running 
time is 0(n\og* s + n log n/ logs) when s working variables are available (and s € f2(logn) n 
O(2 log ™/ log ")). Whenever only O(logn) variables are available, Raman and Ramnath [TS] gave 
an algorithm whose running time is 0(nlog 2 n). More recently Chan [3] provided several lower 
bounds for performing selection with few variables and modified the randomized algorithm of 
Munro and Raman so that it works for any input array. 

There has also been an interest in finding algorithms for geometric problems that use a constant 
number of variables: Given a set of n points, the well-known gift-wrapping algorithm (also known 
as Jarvis march |16j ) can be used to report the points on the convex hull in 0(nh) time using a 
constant number of variables, where h is the number of vertices on the convex hull. Asano and 
Rote S\ and afterward Asano et al. [2] gave efficient methods for computing well-known geometric 
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structures, such as the Dclaunay triangulation, Voronoi diagram and minimum spanning tree using 
a constant number of variables (in 0(n 2 ), 0(n 2 ) and 0(n 3 ) time, respectively) Observe that, since 
these structures have linear size, they are not stored but reported. 

To the best of our knowledge, there is no known method that computes the visibility polygon 
of a given point using few variables. The question of finding a constant workspace algorithm to 
compute the visibility polygon in a polygon of size n was explicitly posed as an open problem 
by Asano et al. PQ. A natural approach to the problem would be to transform the linear time- 
algorithm of Joe and Simpson jTT] to a constant workspace algorithm. However, their algorithm 
cannot be directly adapted to our setting, since their method uses a stack which could contain up 
to fi(n) vertices. 

Results In this paper, we consider a completely different approach. It is easy to realize that the 
differences between V and the visibility polygon of a given point depends on the number of reflex 
vertices of V . For example, if V is a convex polygon, the visibility polygon will be equal to V , but 
we cannot expect this fact to hold when the number of reflex vertices grows. Therefore whenever 
possible we will express the running time of our algorithms not only in terms of n, the complexity 
of V, but also in terms of r and f (the number of reflex vertices in the input and the number of 
them also present in the output, respectively). This approach continues a line of research relating 
the combinatorial and computational properties of polygons to the number of their reflex vertices. 
We refer the reader to 0] and references found therein for a deep review of existing similar results. 

In Section [3] we give an output-sensitive algorithm that reports the vertices of the visibility 
polygon in 0(nf) time using O(l) variables. In Section [4] we show that if we are allowed to use 
0(log r) variables, the problem can be solved in 0(n log r) randomized expected time or 0(n log 2 r) 
deterministic time. 

2 Preliminaries 

Model definition and considerations on input / output precision We use a slight variation 
of the constant workspace model (sometimes also referred as log space), introduced by Asano [3|. In 
this model, an algorithm can use a constant number of variables, and assumes that each variable or 
pointer contains a data word of O(logrt) bits. Implicit storage consumption required by recursive 
calls is also considered a part of the workspace. 

The input of the problem is a polygon V in a read-only data structure. In the usual constant 
workspace model, we are allowed to perform random access to any of the values of the input 
in constant time. However, in this paper we consider a weaker model in which the only allowed 
operation to the input is obtaining coordinates of the next counterclockwise vertex of a given vertex 
of V . This is the case in which, for example, the vertices of V are given in a list in counterclockwise 
order. 

Many other similar models exist in the literature. We note that in some of them (like the 
streaming |10j or the multi-pass model [6]) the values of the input can only be read once or a 
fixed number of times. We follow the model of constant workspace and allow scanning the input 
as many times as necessary. 

The algorithm given in Section [3] uses a weaker model than the constant workspace model 
(since random access to the vertices of the input is not used) . The methods given in Section [1] 
use a logarithmic number of variables, providing another compromise between the space and time 
complexity of this problem. 

Let n be the number of vertices of V. We do not make any assumptions on whether the input 
coordinates are rational or real numbers (in some implicit form). The only operations that we 
perform on the input are determining whether a given point is above, below or on a given line and 
determining the intersection point between two lines. In both cases, the line is defined as passing 
through two points of the input, hence both operations can be expressed as finding the root of 
linear equations whose coefficients are values of the input. We assume that these two operations 
can be done in constant time. Observe that if the coordinates of the input are algebraic values, we 
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Fig. 1. Left: general setting, vertices that are reflex with respect to q are shown with a white point (black 
otherwise). Right: the visibility polygon Vis-p(g) of q; point p is the shadow of vertex V4. 



can express the coordinates of the output implicitly as "the intersection point of the line passing 
through points Vi and Vj and the line passing through Vk and v" (where Vi,Vj,Vk and vi are 
vertices of the input). 

Other definitions The boundary of V is denoted by dP. We regard V as a closed subset of 
the plane. We are also given a point q inside V, from where the visibility polygon needs to be 
computed. The segment connecting points p and q is denoted by pq. We say that a point p £ V is 
visible (with respect to q) if and only if pq C V (otherwise we say that p is not visible) . The set of 
points that are visible from q is denoted by Vis-p(q) and is called the visibility polygon of q. It is 
easy to see that the visibility polygon of q is a closed polygon whose vertices are either vertices of 
V or the intersection point between a segment of V and a ray emanating from q (see Figure [I]) . 

From now on, for simplicity in the explanation, we assume that no line passes through q and 
two vertices of V . We note, however, that our algorithms can be adapted to work in the general 
case with only minor changes. We also need to define what a reflex vertex is in our context. Given 
any vertex Vk, the line £k passing through and q splits V\£k into disjoint components. We say 
that vertex Vk is reflex with respect to q if the angle at the vertex interior to V is more than 7r 
and the vertices Vk-i and Vk+i lie on the same connected component of W. 2 \ 1^ (see Figure [TJ. 
Observe that any vertex that is reflex with respect to q is a reflex vertex (in the usual sense) , but 
the converse is not true. 

Since the point q is fixed, from now on we omit the "with respect to g" term and simply refer 
to these points as reflex. Intuitively speaking, reflex vertices are the vertices where important 
combinatorial changes might occur in the visibility polygon. That is, the positions where the 
polygon boundary can change between visible and invisible. Let r be the number of reflex vertices 
of V . We also define f as the number of reflex vertices of V that are present in Vis-p(<7). Observe 
that we always have r < r < n. 

Given any two points P\,P2 € dV, there is a unique path that travels counterclockwise along 
the boundary of V from p\ to p2- Let Chain(pi,p2) be the set of points traversed in this path, 
including both p\ and P2 (this set is called the chain between p\ and p^)- We say that a chain 
is visible if all the points of the chain are visible from q. A visible chain C = Chain(pi,p2) is 
CCW-maximal if it is visible and no other visible chain starting at p\ strictly contains C. We will 
use the term Chain(p,p) to denote the whole polygon boundary. 

In this paper, we will often use a ray shooting-like operation, which we call RayShooting(p). 
This is a basic operation that given a point p G dV, considers the ray p(q,p) and reports the first 
point of V hit by the ray. This operation can be performed in linear time using constant working 
space, by scanning the edges of V one by one and reporting the point closest to q that intersects 
the ray. Another similar operation that our algorithms will use is computing the shadow of a given 
reflex vertex v. For that case, we will use the FindShadow(p) operation, which gives the shadow 
of p if p is a reflex vertex (p otherwise). 



4 



Luis Barba, Matias Korman, Stefan Langerman, and Rodrigo I. Silveira 



Given a point p € P, let 9(p) be the angle < 9(p) < 2ir that the ray emanating from q and 
passing through p makes with the positive x-axis. Let Vq be the closest point on dP to q, such 
that 9(vq) — 0. It is easy to see that vo is visible and can be computed in linear time. Without 
loss of generality, we will treat vq as a vertex (even though it does not need to be one). Moreover, 
we will assume that the vertices are numbered such that v is on edge v n V\. 

A point p on dP is the shadow of a reflex vertex v if p is collinear with q and v, and p is visible 
from q. Due to the general position assumption, v must be unique and p must be an interior point 
of an edge. As a result, each visible reflex vertex is uniquely associated to a shadow point (and 
vice versa). 

3 An 0(nf) algorithm using 0(1) variables 

In this section we present an output-sensitive algorithm that spends linear time for each reflex 
vertex of the visibility polygon. The idea of the algorithm is to compute maximal visible chains 
as they appear on the boundary of P. Each iteration of the algorithm starts from a point that is 
known to be visible, and finds the last point of the CCW-maximal chain that starts at that point. 
This is repeated until the initial vertex is found again. 

The following lemma characterizes the endpoints of CCW-maximal chains. 

Lemma 1. Let p G dP be a point visible from q, and let C = Chain(p,p') be a CCW-maximal 
chain. Then p' is either (i) equal top, (ii) a reflex vertex ofV, or (Hi) the shadow of some reflex 
vertex of V '. 

Proof. Clearly, all the boundary of P is visible if and only if p' = p. If p =/= p' ', then when going 
on dP from p to p' , p' is the last point found that is visible (i.e. making the chain longer by any 
positive amount would make it include at least one point that is not visible) This can happen for 
only two reasons: either p' is a reflex vertex, or there is some vertex v such that q, v,p' are collinear 
and appear in that order when viewed from q. In other words, if p' is a shadow. □ 

Based on the previous lemma, the algorithm will start from vq and walk on dP in counter- 
clockwise direction, until the endpoint of the current CCW-maximal chain is identified. 

Let p be a visible point and let v r be the first reflex vertex found on dP when going counter- 
clockwise from p. Chain(p, v r ) and the segments qp and qv r define a region that we will denote by 
1Z(p,v r ). The following observation is crucial for the algorithm. 

Lemma 2. Let p be a visible point and let v r be the first reflex vertex encountered on dP when 
going from p in counterclockwise direction. Then v r is visible if and only iflZ(p,v r ) contains no 
vertex from P , or equivalently, if and only if no edge ofV crosses the segment qv r . 

Proof. Vertex v r is visible if and only if there is no edge of P crossing qv r . Such edge exists if only 
if there is at least one vertex inside lZ(p,v r ). This is because p is visible, thus qp is not crossed 
by any edge, and P is a simple polygon, thus Chain(p, v r ) cannot be crossed by an edge either. 
Therefore the only way for an edge to enter lZ(p,v r ) is through the line segment qv r . The result 
follows. □ 

Let FindNextReflexVertex(p) be a routine that returns the first reflex vertex v r found on 
dV counterclockwise from p in 0(n) time. If no such vertex exists the remaining portion of V is 
convex (and thus all vertices until Vq can be reported as visible). Otherwise, vertex v r is identified 
and we need to find the endpoint p' of the chain. To that end we will determine whether lZ(p, v r ) 
is empty, or equivalently, whether RayShooting(iv) = v r . 

If so, we have p' = v r and we can report all vertices in the chain Chain(p, v r ). Then the next 
visible chain must start at the shadow of v r , which can be found with the FindShadow(iv) 
operation. Otherwise, from Lemma[l]we know that v r is not visible and p' must be the shadow of 
some reflex vertex in Chain(w r ., vq). To find that reflex vertex we need one more observation. 
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Fig. 2. When the next reflex vertex, v r in the figure, is not visible, the end of the current visibility chain 
(p') can be found by walking on dV, in order to find the vertex with smallest angle (with respect to qp) 
that lies in lZ(p, v r ) (shaded region); in the example, that vertex is tij. 



Observation 1 Let p be a visible point and let v r be the first reflex vertex encountered on dV 
when going from p in counterclockwise direction. If v r is not visible, then the CCW -maximal chain 
starting at p ends at the shadow of the reflex vertex in lZ(p, v r ) with smallest CCW-angle with 
respect to the segment qp. 

Therefore the algorithm only needs to find the reflex vertex with smallest angle that lies inside 
1Z(p,v r ), which can be done in 0{n) time by simply walking on &P and keeping track of the 
intersections with the line segment qv r . Figure [2] illustrates a step of the algorithm. 

Observe that in all cases, the total time spent in finding the end of the current chain and the 
beginning of the next one is 0(n). Hence, it follows that the total time of the algorithm is now 
0(nf), where f is the number of vertices reflex w.r.t. q of V that are part of the visibility polygon 
of q. 

Moreover, this algorithm can be implemented such that only one pass through the polygon 
vertices is needed for each reflex vertex that appears in the output (plus one initial pass for finding 
Vrj). Algorithm [l] sketches how the algorithm would look. Note also that such implementation 
requires only 3 variables and 1 bit of additional space. 

Theorem 1. The visibility polygon of a point q can be computed in 0(nf) time using constant 
workspace, where f is the number of reflex vertices that are part of the visibility polygon. More 
precisely, it can be computed in (f + 1) passes through the input. 

Proof. It is easy to see that the above algorithm computes the visibility polygon of q, reads the 
input at most f+1 times, and uses a constant number of variables. Indeed, if one avoids calling 
procedures, and reuses variables once they are not needed anymore, the above algorithm can be 
implemented with only three variables. 

Regarding the running time, other than reporting vertices that are in the visibility polygon 
of q, the algorithm only spends extra time in finding the next visible reflex vertex, which takes 
O(n) time. If it is not visible, then the last point of the current visible chain must be found, which 
takes also 0(n) time. Such point is either a visible reflex vertex or its shadow. In either case, the 
O(n) time spent in each loop is charged to the reported reflex vertex. Finally, observe that a reflex 
vertex cannot be charged twice because the algorithm always proceeds by increasing angle with 
respect to the abscissa. □ 



4 An O(nlogr) algorithm using O(logr) variables 
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Algorithm 1 Output-sensitive algorithm for visibility polygon 



p vq, r < hoo 

repeat 

Walk starting from p until next reflex vertex before Do 
r index of reflex vertex found (or +00 if none found) 
if r < +00 then 

Walk from v r until p, keep track of vertex with smallest angle inside lZ(p, v r ) 
j <s— vertex with smallest angle in lZ(p, v r ) (j—+oo if none found) 
if j — +00 then 
(* v r is visible *) 

Walk and report all points from p until finding v r , including v r (but not p) 
Let v s be the shadow of v r 
p v s 
else 

(* Vj found, v r is not visible *) 
Let p' be the shadow of Vj 

Walk and report all points from p until finding p' , including p' (but not p), and Vj 

p <- Vj 

end if 
end if 
until r = +00 

Report all vertices between p and vq (but not p); report vo only if it is a vertex of V 



c 
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In this section we take a completely different 
approach to solve the problem. We consider the 
visibility problem in polygonal chains (instead 
of the whole polygon) and use a divide and 
conquer approach to split the problem into two 
subproblems. For this purpose, we must adjust 
the visibility definitions to make them consis- 
tent for chains. 

Let C = {p , . . . ,Pi, . . . ,p m } be a polyg- 
onal subchain of the boundary of V, such 
that po,p m are both visible points on the 
boundary of V and Pi is a vertex of V, 1 < 
i < 771—1. Assume without loss of general- 
ity that a — 0(po) < 9(pm) = and let 
"Pe = {q,Po,Pi,--- ,Pm, q} be the simple poly- 
gon contained in V obtained by joining q with 
both endpoints of C; see Figure [3] We say that 
a point 1 on a polygonal chain C is C-visible 

(from q), if the segment qx is completely contained in the polygon Vc- Let Vise(<z) 
x is C-visible} be the set of all C-visible points of C. 



PO 



(Ik'' 



Fig. 3. Polygonal chain C and its associated polygon 
Vc- Point x is C-visible, y is not. 



{ieC 



Lemma 3. Let C = {po, ■ ■ ■ ,p m } be a polygonal chain such that Po,p m a-re both visible points 
on the boundary of V , and let x be a visible point inside C lying on the edge pjPj+\ (for some 
j < m). If d = {po, . . . ,pj,x},C 2 = {x,p j+ i, . . . ,p m }, then Vis c (q) = Vis Cl (?) U Vis C2 (g) and 
Vis Cl (q) n Vis C2 (q) = x. 

Proof. We will first show that Vise (q) C Visd (q) U Visc 2 (<?) ■ Let p be a C-visible point belonging to 
C\. By definition, the segment qp is contained in the polygon Vc — {q,Pa, ■ ■ ■ ,Prm <?}• In particular, 
no point of the subchain C\ crosses the edge qp, hence it is Ci-visible. Analogously, we obtain that 
points that are C-visible and belong to C2 are C2-visible. Finally, observe that a point can only be 
visible in both subchains if it belongs to both of them (hence p G Visci (q) H Visc 2 (q) , implying 
that p = x). 
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To prove the other inclusion, let v G Visc 1 (?) such that x ^ v. Assume that qv intersects 
Ci at a point w, since otherwise v would be C-visible. Recall that qv is completely contained in 
the polygon Vd = {q,Po, ■ ■ ■ ,Pj,x,q}, hence w lies on the interior of Vc x - Observe that p m does 
not belong to Vd, hence by the Jordan curve Theorem, when we walk on C from p m to w, we 
will either cut C\, or one of the segments qpo or qx. Whichever situation occurs, we will obtain 
a contradiction: in the first case we will find that the polygon is not simple, and in the second 
case we will reach a contradiction with the visibility properties of either pq or x. The proof when 
v e Visc 2 (<7) is analogous. □ 

Using a similar argument it is easy to see that if x is C-visible, and both Po,p m are visible from 
g, then a; is a visible point of V. Thus Lemma [3] allows us to divide the problem of finding Vis-p(g) 
into a series of subproblems that can be solved independently. 

Let A(C) be the interior of the cone with apex q defined by the rays going from q and passing 
through both endpoints of C (and crossing the interior of the chain). Our divide and conquer 
algorithm is as follows. On each step of the algorithm we choose a random reflex vertex z inside 
the cone A(C), we then perform a ray shooting query to find the first point x on C in the direction 
of z. We split the polygonal chain C at x, thereby obtaining two subchains C\,Ci. We repeat the 
process recursively first on C\ and then on C2, until C is split into a series of subchains, each 
containing no reflex vertex inside their associated cone. Since the changes in visibility occur only 
at reflex vertices, the visible vertices inside each split subchain can be reported as visible without 
any additional information from the other subchains, as long as they are reported in the order in 
which they are found; see Figure [4] 




Fig. 4. Left: Split of C into two subchains C\,Ci using a visible point x in the direction of a random reflex 
vertex z. Right: A polygonal chain C with no reflex vertices inside the cone A(C), note that only one 
subchain of C is visible. 



The main algorithm presented in this section is summarized in Algorithm [2] The subroutine 
Report VisibleChain(cc) takes a visible point x on C, walks from x toward p m and reports every 
vertex until finding an edge e of C intersecting the ray p(q,p m ), once found that intersection is 
reported and the subroutine ends. 

Lemma 4. Algorithm^ reports every visible edge of Vise (q) in counterclockwise order. 

Proof. Let C = {po, . . . ,p m } such that po,p m are both visible points of V . The proof goes by 
induction on k, the number of reflex vertices inside the cone A(C). If k = then only one subchain 
of C is visible, and finding it only requires a ray shooting query to find a point x being the shadow 
of pq , followed by a walk from x toward p m on the path to report all visible points inside the cone. 
Since the walk is performed in counterclockwise order, the output is also reported in the same 
order. 

If k > 0, then the algorithm creates two subchains Ci, C2 of C, such that the cones A(C\), Mf^-i) 
contain at most k — 1 reflex vertices each. By the induction hypothesis the output of the recursion 
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Algorithm 2 Given a polygonal chain C = {po, . . . ,p m } such that Po,p m are both visible points 
of V, algorithm to compute Vis c (q) 

k number of reflex vertices of C inside the cone A(C) 
if k = then 

s •(— FindShadow(po) 
if s 7^ po then 

Report po 
end if 

ReportVisibleChain(s) 
else 

r <— random number in {1, . . . , k} 
z the r-th reflex vertex of C inside A(C) 
x RayShooting(z) 
Call Algorithm [5] on Ci = {po, ■ ■ ■ ,x} 
Call Algorithm on C2 = {x, ... ,p m } 
end if 



on C\,C2 is exactly Vis^ (q), Visc 2 (<?) respectively, each reported in counterclockwise order. Thus 
by Lemma [3] it follows that Algorithm [2] reports the set Vise (q) correctly, with vertices given in 
counterclockwise order. □ 

Lemma 5. The expected running time of Algorithm^ is 0(n log r). 

Proof. The running time of Algorithm [2] can be analyzed in a similar way to the one of quicksort. 
On each step, a random pivot is chosen and it is used to split the chain C, such that each subchain 
obtained contains in expectation a fraction of the reflex vertices. Thus if we look at the recursion 
tree, it is easy to see that on each level 0(n) operations are needed and the expected depth of the 
tree is O(logr), resulting in an expected running time of 0(n log r). 

The exact details of the proof are as follows: the recursion of the algorithm depends on two 
variables n, k, the number of vertices on the polygonal chain C — {po, . . . ,p m }, and the number of 
reflex vertices inside the cone A(C), respectively. Imagine that the reflex vertices inside A(C) are 
in a list L, sorted in counterclockwise order according to their angles with respect to qpo. Let Xj 
be an indicator random variable such that 



X r = 



1 If the z-th reflex vertex on L was chosen by the algorithm 
Otherwise 



Since we chose the pivot vertex at random, then Pr[Xi = 1] = 1/k. The running time of our 
algorithm can be expressed as follows . 

k 

T(n, k) = ^r,Xi (T(hi,i) +T(n-h h k- i)) + 0{n). 

i=l 

Where hi is the number of vertices on the chain between vq and the i-th reflex vertex of L. We 
now analyze the expected running time, therefore 

k 

E[T(n, k)} = E[X t (T(h t , i) + T(n -h u k- *))] + 0(n) = 
»=i 

1 k 

- V E[T(hi, i)] + E[T(n -K,k-i)} + 0{n) 

i=l 

We claim E[T(n, k)] € 0{n log r) and we prove it by substitution, suppose E[T(to,j)] < am log j 
for all to < n, j < k, thus 

k 

E[T(n, k)} < - Y hi logi + (n - hi) log(fc - i) 

i=l 
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fc k/2 

E[T(n, fc)] < - V n log(max{«, fc - i}) + 0(n) = — V logffc - i) + 0(n) 
fc ^— ' fc ' 

i=l i=l 

Using the integral bound we obtain that 

E[T(n, fc)] < -y f - log -J + 0(n) = an log - + 0(n) = 

anlogfc — (an log 2 + 0(n)) . 

Thus by picking a sufficiently large to dominate the constant on the 0(n), E[T(n, fc)] < anlogfc 
and since fc is the number of reflex vertices inside the cone, we have fc < r. We conclude that the 
expected running time of Algorithm [2] is 0(n log r). □ 

Theorem 2. T/ie visibility polygon of point q can be computed in 0(n log r) expected time using 
O(logr) variables. 

Proof. In each step of the recursion a constant number of variables are needed. Hence, the number 
of words used by the algorithm is proportional to the depth of the recursion. In order to avoid 
using an excessive amount of memory, we use a standard trick of restarting the algorithm whenever 
the recursion tree becomes too deep, as follows. 

We introduce a small variant of the algorithm, to be able to verify that the algorithm never 
uses more that O(logr) words. It is easy to see that the recursion tree of the Algorithm [2] is a 
random binary tree, therefore its expected depth is clogr, where c s» 4.311 [7]. Let A denote the 
depth of the recursion tree of Algorithm [2] on an arbitrary run and recall that from Markov's 
inequality we have Pr[X > 2c log r] < 1/2. Let wbea constant denoting the number of words used 
by the algorithm on each step of the recursion. 

If at some point in the execution of the algorithm, the program wants to allocate more that 
2wclogr words, we stop and start again from scratch. From Markov's inequality, the expected 
number of times we need to restart the algorithm before it finishes without exceeding the 2u;c log r 
limit is one. Since the expected number restarts is bounded by a constant, the (asymptotically 
speaking) expected running time of the algorithm is unchanged. 

Finally note that by restarting the algorithm we might report many times some vertices of 
Vis-p(g). In order to avoid so, at any point in the execution of the algorithm, we store the last 
vertex v that we reported (this value is the only one that is not reset when the algorithm is 
restarted). Whenever the algorithm wants to report a point p, we compare the angles 9(p) and 
6(v) and only report p if 9(p) > 6(v). Since the vertices are reported in counterclockwise order, 
then no point is reported twice. □ 

We now describe a deterministic variant of Algorithm [5] that runs in 0(n log 2 r) time using 
O(logr) words. Let 7Z — {6(vi) : Vi is a reflex vertex inside A(C)} and let n, fc be the number of 
vertices of C and the cardinality of 1Z, respectively. Recall that in steps [9] and 10 of Algorithm [2] 



we are choosing a random reflex vertex v inside A(C) and using it to split the chain C into two 
subchains Ci,C2, such that the number of reflex vertices inside A(C\) and A(C%) is balanced in 
expectation. The next variant of Algorithm [2] replaces the random selection with a deterministic 
selection algorithm that guarantees a balanced split, albeit at a slight increase in the running 
time. The algorithm is based on the approximate median pair algorithm proposed by Raman and 
Ramnath [15] . which allows to compute an approximation of the median of a set of m elements 
using O(logm) variables in O(logm) passes. 

Theorem 3. The visibility polygon of point q can be computed in 0(n log 2 r) time using O(logr) 
variables. 

Proof. Let 71, n and fc be defined as before. We describe a variant of Algorithm [2] that replaces the 
random selection with a deterministic selection algorithm that guarantees a balanced splitting. 

Let z,z' be two elements of 1Z such that z < z' , and let Isr be the list obtained by sorting 
TZ. We say that z, z' is an approximate median pair if z and z' partition L-ji in three parts, each 
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of size at most k/2. Note that if z, z' is an approximate median pair, then at least one of z or z' 
splits L-ji into two sets each of size at most 2k /3. Moreover, deciding which one of them to use 
requires a single extra pass on TZ. Thus we say that every approximate median pair induces an 
even split of TZ. 

We use the selection algorithm by Raman and Ramnath |15j to compute an approximate 
median pair of a set of m elements using O(logm) variables in O(logm) passes. Since \TZ\ — r and 
the cost of each pass is 0(n), we can implement an algorithm to find an approximate median pair 
of TZ in 0(n\ogr) time using O(logr) variables. Using the even split induced by this algorithm, 
we can find a reflex vertex inside A(C) that splits the chain C into two subchains, each containing 
at least fc/3 reflex vertices inside their associated cone. Thus the recursion tree has always depth 
O(logr). Moreover, each level of the algorithm costs O(nlogr), therefore the result follows. □ 

5 Closing Remarks 

One expects that the running time of a constant workspace algorithm increases when compared 
to other models in which memory is not an issue. In most geometric problems, the increase in the 
running time is almost equal to the reduction in memory space 2 3 . The algorithm given in Section 
[3] follows a similar pattern, since in the worst case the time-space product matches the one of Joe 
and Simpson (from 0(n 2 ) to 0(nf)). However, notice that the improvement of the algorithm 
given in Section [4] is much better, since the time-space product becomes 0(n log n) x O(logr) = 
O(nlognlogr). 

If we compare the two methods given in this paper, we obtain a linear reduction in the running 
time by increasing the memory space by a logarithmic factor. This is due to the fact that this space 
allowed us to use more powerful techniques, such as divide and conquer. It would be interesting 
to study if the same result holds for other geometric problems (such as computing the Voronoi 
diagram, Delaunay triangulation or convex hull). 
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